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A CHANDRA STUDY OF THE CIRCINUS GALAXY POINT-SOURCE POPULATION 
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ABSTRACT 

We have used the Chandra X-ray Observatory to resolve spatially and spectrally the X-ray emission from the 
Circinus Galaxy. We report here on the nature of the X-ray emission from the off-nuclear point sources associated 
with the disk of Circinus, which make up w 34% of the total 0.5-10 keV emission. We find that many of the 
serendipitous X-ray sources are concentrated along the optical disk of the galaxy, although few have optical 
counterparts within 1" of their X-ray positions down to limiting magnitudes of my = 23-25. At the distance of 
Circinus (w 3.8 Mpc), their intrinsic 0.5-10 keV luminosities range from w 2 x 10"'^ erg s~^ to « 4 x 10"''^ erg s~' . 
One-fourth of the sources are variable over the duration of the 67 ks observation, and spectral fitting of these 
off-nuclear sources shows a diverse range of spectral properties. The overall characteristics of the point sources 
suggest that most are X-ray binaries and/or ultra-luminous supernova remnants within Circinus. 

We are able to analyze the two strongest off-nuclear sources in greater detail and find both to have remarkable 
properties. The average X-ray luminosities of the two sources are 3.7 x lO"*^ erg s"' and 3.4 x lO-'^ erg s"'. The 
former displays large and periodic flux variations every 7.5 hr and is well fit by a multicolor blackbody accretion- 
disk model with ?]„ = 1.35 keV, properties consistent with an eclipsing > 50 Mq black-hole binary. The latter 
appears to be a young supernova remnant, as it coincides with a non-thermal radio counterpart and an Ha-detected 
H II region. This source exhibits both long-term (« 4 yr) X-ray variability and a 6.67-6.97 keV iron emission-line 
blend with a 1 .6 ke V equivalent width. These two objects further support the notion that super-Eddington X-ray 
sources in nearby galaxies can be explained by a mixture of intermediate-mass black holes in X-ray binaries and 
young supernova remnants. 

Subject headings: galaxies: active — galaxies: individual (Circinus) — X-rays: galaxies - X-rays: binaries 



1. INTRODUCTION 

The study of X-ray emission from star-forming galaxies has 
long lagged behind studies in more traditional wavebands (e.g., 
ultraviolet, optical and radio) because of sensitivity and reso- 
lution limitations. These constraints largely restricted detailed 
investigations to sources in the Local Group, where only small 
pockets of ongoing star formation exist. The Einstein X-ray 
Observatory was the first to resolve the X-ray source popula- 
tions in nearby galaxies (e.g., Helfand 1984; Fabbiano 1989), 
opening the study of luminous X-ray binaries (XBs) and su- 
pernova remnants (SNRs) to sources outside our own Galaxy 
(and the Magellanic Clouds) where there is generally less am- 
biguity in the distance to the source and hence its intrinsic lu- 
minosity. The brightest point sources in star-forming galaxies 
were found to be significantly more luminous on average than 
the brightest sources in the Milky Way (e.g., Fabbiano 1989). 
More sensitive X-ray observations of normal and star-forming 
galaxies with ROSAT and ASCA have shown that variable, off- 
nuclear point sources with super-Eddington X-ray luminosi- 
ties (i.e., Lx > 2 X lO"*^ erg s"' for M = I.4M0) are common 
but generally not well understood (e.g., Colbert & Mushotzky 
1999; Makishima et al. 2000; Roberts & Wai'wick 2000; Lira, 
Lawrence, & Johnson 2000). 

The large improvements in spatial resolution and sensitivity 
recently afforded by Chandra, however, have begun to change 



this picture. Chandra not only offers detailed studies of the 
brightest point sources in neighboring galaxies, but also studies 
of these sources over a broad range of galaxy types and lumi- 
nosity classes, providing large samples of X-ray point sources 
that can be related by luminosity and spatial distribution to 
other properties of these galaxies such as gas mass, star for- 
mation rate and morphology (e.g., proximity to spiral arms, 
H II regions or globular clusters). Recent Chandra observa- 
tions of M31 (Garcia et al. 2000), M81 (Tennant et al. 2001), 
M82 (Griffiths et al. 2000), NGC 3256 (R Lira et al., in prepa- 
ration), NGC 4038/39 (Fabbiano, Zezas, & Murray 2001), and 
NGC 4647 (Sarazin, Irwin, & Bregman 2000) are just a few of 
the diverse examples. 

Here we perform a high-resolution X-ray study of the mas- 
sive spiral known as the Circinus Galaxy (hereafter Circinus). 
Because of its proximity to the Galactic plane (b = -3.°8), Circi- 
nus lay hidden until the 1970s. Freeman et al. (1977) found 
that Circinus lies at a distance of 3.8 ± 0.6 Mpc"* with an in- 
clination of ^65° to the line-of-sight and appears to be lo- 
cated within a Galactic "window" with a visual absorption 
of Ay = 1.5 ±0.2 and a neutral hydrogen column density of 
A^H = (3.0±0.3) X 10^' cm-2 [whereas neighboring regions typ- 
ically have Ay = 3.0 and A^h = (5-10) xlO^' cm"-; Schlegel, 
Finkbeiner, & Davis 1998; Dickey & Lockman 1990].'* Since 
its discovery, numerous lines of evidence have been presented 
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indicating that this galaxy hosts a Seyfert nucleus: a prominent 
[O III] ionization cone (Marconi et al. 1994; Wilson et al. 2000), 
near-IR polarized bipolar scattering cones (Ruiz et al. 2000), a 
coronal line region (Oliva et al. 1994; Maiolino et al. 2000), po- 
larized broad Ho; (Oliva et al. 1998; Alexander et al. 2000) and 
strong X-ray iron Ka Unes (Matt et al. 1996; Sambruna et al. 
2001b). Additionally, this galaxy displays a moderate level of 
star formation, as indicated by Ha imaging of a complex net- 
work of H II regions (e.g., ~ 10^"^ yr old star-formation rings 
at ^ 40 and ^ 200 pc from the nucleus; Marconi et al. 1994; 
Wilson et al. 2000). Past X-ray studies have concentrated on 
the obscured nuclear source, with the only spatially resolved 
pTS-Chandra observation of Circinus being a 4 ks ROSATHigh- 
Resolution Imager (HRl) detection of the nucleus and one off- 
nuclear point source (Guainazzi et al. 1999). 

In this paper, we study the properties of 16 X-ray 
point sources coincident with Circinus which have fluxes 
■Pb.s-io.o kev > 7 X 10"'^ erg cm"^ s"'. Details of the observation 
and reduction procedure are outlined in §2, and source identifi- 
cations are discussed in §3. In §4 and §5 we discuss timing and 
spectroscopic analyses of the X-ray sources. We comment on 
the luminosity function of the point sources in Circinus in §6, 
and, finaUy, we discuss the implications of our findings in §7. 

2. OBSERVATIONS AND DATA REDUCTION 

Circinus was observed with the High-Energy Transmission 
Grating Spectrometer (HETGS; C. R. Canizares et al. in prepa- 
ration) on 2000 June 6-7 with ACIS-S (G. R Garmire et al., in 
preparation) in the focal plane. The HETGS provided both a 
zeroth-order image of Circinus on the ACTS S3 CCD and high- 
resolution dispersed spectra from two transmission gratings, the 
Medium-Energy Grating (MEG) and the High-Energy Grating 
(HEG), along the ACIS S0-S5 CCDs. The Circinus nucleus 
was placed near the aim point on the S3 chip and was observed 
using a subarray window of 600 rows to reduce pileup in the 
zeroth-order image.'' This particular configuration provided a 
frame time of 2.1 s and reduced pileup in the nucleus (i.e., the 
brightest X-ray source in the field) to « 3% and in the off- 
nuclear sources to < 1% (see Sambruna et al. 200 1 a). ^ Events 
were telemetered in Faint mode, and the CCD temperature was 
-120° C. Chandra is known to experience periods of especially 
high background which are most pronounced on the S3 chip,** 
but such "flares" are much less frequent when the HETGS is in- 
serted.' We found that no background "flares" occurred during 
this observation. 

Analysis was performed on reprocessed Chandra data (2000 
December reprocessing) primarily using the CI AO V2.0 soft- 
ware provided by the Chandra X-ray Center (CXC), but also 
with FTOOLS and custom software. We removed the 0"5 pixel 
randomization, performed standard ASCA grade selection, and 
excluded bad pixels and columns. The observation was split 
into two continuous segments separated by 1,270 s because the 
HETGS was not fully inserted during the first segment (^ 0.°3 
from full insertion). While the HETGS misalignment had a pro- 

* Pileup is the coincidence of two or more photons in a pixel per ACIS frame that are counted by the instrument as a single, higher energy event. The presence of 
pileup results in photometric inaccuracy as well as in distortions to the HRMA+ACIS point spread function and to ACIS spectra. 

^ The frame time is the fundamental unit of exposure for readout of the ACIS CCDs. The default frame time is 3.2 s, but it can range from 0.2-10.0 s, depending on 
the number of ACIS chips and subarray used. 

* See http://asc.harvard.edu/cal/Links/Acis/acis/Cal_prods/bkgmd/10_20/bg201000.html. 
' See http://space.mit.edu/HETG/flight/status.html. 

See http://space.mit.edu/HETG/technotes/offset_effect_000619.txt. 

We note that there may be another diffuse emission component arising from imresolved point soitfces and/or ionized gas associated with the optical star-forming 
disk. However, this X-ray emission is diluted by the bright reflection halo that extends out to large radii (> 1' from the nucleus). 



found effect on the gratings data, changes to the zeroth-order 
imaging and spectroscopy should be negligible."' The total 
net exposure times for these two segments were 6,933 s and 
60, 153 s, respectively. The two datasets were merged (67,086 s) 
to enhance source detection and timing analysis. Since the tele- 
scope did not move between the two exposures, we also ex- 
tracted spectra using the fuU exposure. Event PI values and 
photon energies were determined using the latest gain files ap- 
propriate for the observation. The X-ray spectra were analyzed 
using XSPEC (Arnaud 1996). Unless stated otherwise, spectral 
parameter errors are for 90% confidence assuming one param- 
eter of interest. 

The X-ray fluxes and absorption-corrected luminosities for 
all sources with more than 90 counts were calculated from 
spectral fitting using XSPEC. For sources below this count 
limit, fluxes and luminosities were estimated assuming an av- 
erage power-law spectrum with photon index F = 1.76 and 
A^H = 1.5 X 10^^ cm"^ as determined from the best-fit values 
to sources above the 90 count limit. These flux and luminos- 
ity estimates are assumed throughout the paper unless stated 
otherwise. The scatter in the observed values of V and A^h sug- 
gests that there is likely to be some systematic uncertainty in 
our low-count flux and luminosity estimates. To assess this 
potential uncertainty, we replaced the average F and A^h val- 
ues with the entire range of values for sources with more than 
90 counts, rather than the average values. This led to devia- 
tions of no more than 36% above and 24% below the average 
flux and 54% above and 15% below the average absorption- 
corrected luminosity. We also note that the average value of A^h 
used above is much larger than the Galactic absorption column 
(3 X 10^' cm"^) and has a large dispersion (8 x 10^' cm"^), im- 
plying that there is likely to be a wide range in A^h values among 
individual sources. 



3. SOURCE DETECTION AND DISTRIBUTION 

Figure 1 shows an HETGS zeroth-order "false color" image 
of Circinus. The colors red, green, and blue represent the soft 
(0.5-2.0 keV), medium (2.0^.0 keV), and hard (4.0-10.0 keV) 
X-ray emission, respectively. Prior to combination, each color 
image was smoothed with an adaptive kernel algorithm (Ebel- 
ing. White, & Rangarajan 2001) which permits the simultane- 
ous viewing of compact and extended sources. The smoothed 
images each have a signal-to-noise ratio of 2.5 per smoothing 
beam and were corrected for exposure. Several features in the 
HETGS image are evident: (1) a central, bright point source, 
(2) a circumnuclear diffuse reflection component, (3) a soft X- 
ray plume coincident with the known ionization cone, and (4) 
several off-nuclear point sources. The point-source distribu- 
tion is elongated in the same direction as the projected major 
axis of the optical disk (a position angle of about -45° East of 
North), with the majority of the sources concentrated towards 
the nucleus. We see a prominent division of the circumnuclear 
reflection halo in Figure 1 into hard (blue) and soft (red) X-ray 
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Fig. 1 . — An HETGS zeroth-order "false color" 80" X 80" image of Circinus with red, green, and blue representing 0.5-2.0 keV, 2.0^.0 keV, and 4.0-10.0 keV 
emission, respectively. Prior to combination, each color image was smoothed with an adaptive kernel algorithm (Ebeling, White, & Rangarajan 2001) to enhance 
diffuse features. The emission features depicted are (1) a central point source, (2) a diffuse circumnuclear reflection halo, (3) a soft X-ray plume coincident with 
the known ionization cone, and (4) several off-nuclear point sources. Most of the off-nuclear sources appear to lie along the projected major axis of the galaxy with 
a position angle of about —45° East of North). The dashed lines represent the opening angle of the [O III] and Ha detected ionization cone (Marconi et al. 1994; 
Wilson et al. 2000). The two brightest off-nuclear X-ray sources, CG X-1 and CO X-2, are marked. 



emission which appears to be perpendicular to the optical disk 
of the galaxy. This division is indicative of intrinsic absorption 
and implies that there is a large gas and dust reservoir lying in 
the plane of the galaxy near the nucleus. This is consistent with 
the large absorption suggested from infrared observations of the 
nuclear region (Ruiz et al. 2000). 

The first three features of Circinus are discussed in Sam- 
bruna et al. (2001a, b). In the following analysis, we focus on 
the zeroth-order imaging and moderate-resolution spectroscopy 
of the fourth component. We used the CI AO source detec- 
tion tool wavdetect (Freeman et al. 2001) to identify all 
point sources above a threshold significance of 10"^. A total 
of 16 sources were found. Table 1 lists these sources along 
with their 0.5-10 keV fluxes, luminosities and spectral char- 
acteristics (see §5). Point-source counts and spectra were ex- 
tracted using the 95% encircled-energy radii at 1.5 keV (4 pix- 
els on-axis and up to 7-8 pixels 2' off-axis)'^ for all sources 
except CXOU J141309.2-652017, CXOU J141310.0-652021, 
CXOU J141310.3-652017, and CXOU J141310.4-652022. 
For these sources, the 90% encircled-energy radii were used 
to limit the amount of contamination from the diffuse circum- 
nuclear emission component (for details on this diffuse compo- 
nent, see Sambruna et al. 2001a). The background-subtracted 
counts measured from the 90% and 95% encircled-energy 
radii were multiplied by 1.11 and 1 .05 to obtain the aperture- 



corrected count values given in column 2 of Table 1. Since 
some sources are embedded in diffuse emission well in excess 
of the typical 0.07 cts pixel"' background for this observation, 
we extracted background counts using circular or annular re- 
gions of 100-1000 pixels adjacent to or surrounding the source 
extraction region. Our absorbed, background-subtracted point- 
source detection limit is 7 counts, corresponding to an absorbed 
flux of ^b.5-iokeV = 7 X 10~'^ erg cm~^ s~' or an absorption- 
corrected luminosity of Lo.s-io keV = 2 x 10^^ erg s~' (assuming 
an average power-law spectrum with photon index F = 1.76 and 
A^H= 1-5 X 10^2 cin-2). 

The optical nature of these serendipitous X-ray sources was 
assessed using HST WFPC2 (F502N, F547M, F606W, F656N 
and F814W) archival images; see Table 2 for details. The im- 
ages were processed through the Space Telescope Science In- 
stitute pipeline.'^ When multiple images in a given filter were 
available, the images were combined to reject cosmic rays using 
the I RAF task crre j. The automatic star-finding algorithm 
daof ind in DAOPHOT (Stetson 1987) was used to generate a 
preliminary list of stars detected in each chip for each filter. We 
then performed aperture photometry to determine each star's in- 
strumental magnitude using 2-pixel radius apertures in the WF 
chips and 5-pixel radius apertures in the PC chip. Aperture cor- 
rections to a 0."5 radius were then determined from bright stars 
on each chip. 



E. D. Feigelson 2001, private communication; see also http://asc.harvard.edu/cal/Links/Acis/acis/Cal_prods/psf/analysis.html 
See http://www.stsci.edu/documents/dhb/web/DHB.html. 



Table 1 

X-RAY Sources in the ACIS-S3 Image of Circinus 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


CXOU J Source 


Counts 


^con 


my 


Fx 


Lx 


log(fx/iv) 


AfH 


r 


Comments 


141257.8-651425 


38.0 ± 10.9 


97.3 


> 20.7 (DSS) 


0.040+ 


0.116+ 


> 0.25 








141303.3-652043 


17.4 ± 6.5 


42.6 


22.1 (HST) 


0.019+ 


0.053+ 


0.49 








141305.5-652031 


21.3 ± 7.0 


38.6 


> 25.3 (HST) 


0.023+ 


0.065+ 


> 1.85 








141307.5-652106 


7.7 ± 3.9 


8.4 


22.2 {HST) 


0.007+ 


0.020+ 


0.10 






Variable (flare) 


141309.2-652017 


99.1 ± 12.0 


43.8 


> 23.4 {HST) 


0.10 


0.24 


> 1.73 


0.99 (< 4.87) 


1 A ^+1 36 

i-44!i:g 




141309.7-652211 


11.9 ± 4.6 


29.8 


> 20.7 (DSS) 


0.011+ 


0.032+ 


> -0.31 








141310.0-652021 


6571.2 ± 81.4 


85.6 


^ 19.0 {HST) 


12.3 


23.3 








Nucleus 


14131U.U— 0DZU44 (CLr A-2) 


1 ^lA C _1_ '3'7 1 




ZZ.l {HSl) 


1 11 
1.33 


1 A1 

3.43 


1 (LI 

z.oi 


1 1/140.29 
1-14-0.26 


1 -77+0.21 
'-0.22 


Long-term variable, 6.9 kcV emission line, 
















emoeciueci in riii region i^naj j 










U.Ul / 


U.U^o 










141310.1-652029 


166.0 ± 15.0 


93.2 


> 25.3 {HST) 


0.17 


0.73 


>2.28 


3.37!i-^» 


2.03!ii| 




141310.3-652017 


342.5 ± 20.6 


85.1 


> 25.3 (i/5D 


0.26 


1.40 


> 1.94 


1.26!g:i 


'T-l+O.Sl 

^- '^-0.52 




141310.4-652022 


792.7 ± 29.7 


2.3 


> 25.3 {HST) 


1.08 


2.09 


>3.08 


0.9lt«i 


0.68!o:2i 


Variable (rise) 


141312.2-652007 


186.0 ± 14.9 


0.9 


20.1 (DSS) 


0.21 


0.52 


0.73 


1.64 « 4.22) 


] O-7+0.80 
i-J '-0.84 


Variable (flare) 


141312.3-652013 (CGX-1) 


1095.2 ± 34.2 


0.0 


> 25.3 {HST) 


0.90 


3.71 


>3.44 


i -^0-0.22 


2.4o:«:i 


Variable (periodic) 


141312.6-652052 


119.2 ± 12.2 


0.4 


> 25.3 (HST) 


0.12 


0.33 


> 1.65 


1.43 «5.26) 


i-63_o.64 


Variable (flare) 


141312.9-651957 


27.7 ±7.1 


16.8 


> 18.1 (DSS) 


0.029+ 


0.084+ 


> -0.93 








r radius aperture 


11760 ± 152 


4.2 




10.0 


17.2 











Note. — Column 1: Source name given as CXOU JHHMMSS.S+DDMMSS. The final row lists the global properties of Circinus determined within a I'aperture radius of the 
nucleus using software provided by A. Vikhlinin (http://asc.harvard.edu/cgi-gen/contributed_software.cgi). Column 2: Background-subtracted, aperture-corrected 0.5-10.0 keV counts 
accumulated over 67 ks. Aperture photometry was performed using 90% and 95% encircled-energy radii for 1.5 keV, and individual background regions were selected adjacent to 
each source as noted in §3. The standard deviation for the source and background counts are computed following the method of Gehrels (1986) and are then combined following the 
"numerical method" described in §1.7.3 of Lyons (1991). Column 3: The probability that the source count rate is consistent with a constant rate, as derived from the Kolomogrov- 
Smimov test (in percent). Column 4: Visual magnitudes or lower limits as determined from HST or DSS images (see §3). Colunm 5: Observed 0.5-10.0 keV fluxes in units of 
10"'^ erg cm"^ s ' from the best-fit models to the ACIS spectra. Fluxes for sources denoted by a f were calculated assuming an average power-law spectrum with A^h = 1-5 x 10^^ 
cm"^ and photon index F = 1.76 as determined from the best-fit values of sources with more than 90 counts. Column 6: Absorption-corrected 0.5-10.0 keV luminosities in units of 
10^' erg s~' from the best-fit models to the ACIS spectra. Luminosities for sources denoted by a f were calculated with the same assumptions as outUned in Column 5. Column 7: Ratio 
of the 0.5-10.0 keV X-ray flux to the V magnitude optical flux. Colunm 8: Neutral hydrogen absorption colunm density in units of 10 cm~^ as determined from the best-fit models to 
the ACIS spectra. Also listed are the 90% confidence errors calculated for one parameter of interest (Ax^ = 2.7). Column 9: Power-law photon index F as determined from the best-fit 
models to the ACIS spectra. Also listed are the 90% confidence errors calculated for one parameter of interest (Ax^ = 2.7). Column 10: Comments on X-ray variability and optical 
identifications. 
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Tabll 2 

IVAL iMACiLS Ol- ClRC'lM S 



Ha Image 



V Image 



(1) 


(2) 


(3) 


(4) 




Filtert 


Exp. Time (s) 


m\im 


Obs. Date 


WFPC2 F502N 


900, 900 


23.3 


1999 Apr. 


10 


WFPC2 F547M 


60 


23.0 


1999 Apr. 


10 


WFPC2 F606W 


200, 400 


25.3 


1996 Aug. 


11 


WFPC2 F656N 


800, 800 


22.7 


1999 Apr. 


10 


WFPC2F814W 


40 


22.8 


1999 Apr. 


10 




-65^ 20' 16' 



10-5'' 14" 13'" 10.0" 

Right Ascension (2000) 




12.5^ 14" 13'" 12.0'' 

Right Ascension (2000) 



t The WFPC2 filter name denotes the central filter wavelength 
in nanometers and the filter width (N=narrow, M=medium, 
W=wide). 

We aligned the HST images to the HipparcosfTycho astro- 
metric reference frame using three moderately bright stars from 
the Tycho catalog (H0g et al. 2000) to provide absolute astrom- 
etry to w 0."4. The F606W image had approximately the same 
aim point as the other WFPC2 images, but it was rotated by 
130° such that most (but not all) of a 2' field centered on Circi- 
nus was imaged by WFPC2. We filled the slight optical cov- 
erage gap using Q'.'l resolution 5100 A and 7000 A continuum 
images taken with the Superb-Seeing Imager on the New Tech- 
nology Telescope (NTT) in 1994 (Figures 1 and 2 of Marconi 
et al. 1994). The Chandra zeroth-order field-of-view is signifi- 
cantly larger than either the HST or NTT images, so the optical 
nature of the few sources farthest from the center of Circinus 
was evaluated using the Digitized Sky Survey (DSS; Lasker et 
al. 1990). Unfortunately, the coarse resolution (FWHM w 3") 
and shallow depth (limiting ms ~ 20) of the DSS image in this 
region of the sky are not nearly as constraining as the HST im- 
ages. The NTT and DSS images were aUgned with the HST 
images using ten bright, isolated point sources in the field. The 
astrometric accuracy of each image is w 1 "0. 

The only obvious optical counterpart to any of the X-ray 
sources is the nucleus of Circinus, which is offset from the X- 
ray nucleus centroid by 0."7. Given the sUght uncertainty in the 
astrometry of the HST images, this value is consistent with the 
typical astrometric accuracy for reprocessed Chandra data (i.e., 
< Q'.'6)}'^ Aligning the images based on this shift, we found 
that three of the off-nuclear X-ray sources are coincident with 
HST optical counterparts to within 0"5: an my = 22.1 coun- 
terpart with CXOU J141303.3-652043, an my = 22.2 counter- 
part with CXOU J141307.5-652106, and an my = 22.8 coun- 
terpart with CXOU J1413 10.0-652044. Johnson V magnitudes 
for these sources were calculated using SYNPHOT assuming a 
spectral slope determined from their wps u - m^^An color. For X- 
ray sources without HST detections, V magnitude limits were 
calculated using SYNPHOT assuming a spectral slope deter- 
mined from the average mpuA — m^5Ai color for sources in the 
field (i.e., mF8i4-ffJF547 « 1-4). Several of the detected X-ray 
sources, however, fall outside of the range of the HST obser- 
vations. Inspection of the NTT 5100 A continuum and DSS 
images reveals an my = 20.1 counterpart at a distance of 0."8 
from another source (CXOU J141312.2-652007). Johnson V 

^'^ See http://asc.harvard.edu/mta/ASPECT/cel_loc/cel_loc.html. 

Johnson B and R magnitudes were derived from the United States Naval Observatory A2.0 Catalog O and E plate magnitudes using using B = 0-Q.ll9{0-E) and 
« = (Evans 1989). 

The typical X-ray-to-optical flux ratios derived by Maccacaro et al. (1988) for individual source types were based on X-ray observations in the 0.3-3.5 keV band 
with low neutral hydrogen column densities. However, all of the Circinus off-nuclear sources exhibit spectral cutoffs below ~ 1 keV (see §5). Therefore, the flux 
ratios listed here have been corrected for the Galactic absorption column and visual extinction given in § 1. The flux ratios have also been converted from 0.3-3.5 keV 
to 0.5-10.0 keV assuming the average power-law spectrum measured for sources in the Circinus field. Large spectral deviations from our adopted model could result 
in log(Fx/fv) differences of up to ~ 0.5-1.0. 



Fig. 2. — HETGS zeroth-order contours of the two brightest off-nuclear 
sources, CG X-2 (left) and CG X-1 (right), overlaid on HST ¥656 WF4 and 
F606 PC chip images, respectively. The lowest contour level is 4 counts, with 
each successive contour increasing by a factor of 4. CG X-2 coincides with an 
my = 22.2 companion embedded in a probable H II region, as indicated by the 
diffuse emission surrounding the optical point source in the HST Ha image. 
CG X-1 has no optical counterpart to my > 25.3. 

magnitudes and lower limits for sources covered by the NTT 
and DSS images were calculated using SYNPHOT assuming 
a spectral slope determined from the average B—R color for 
sources in the field (i.e., B-R k, 1.25).'^ 

In Figure 2, we show contours of the two brightest 
off-nuclear X-ray sources superimposed on HST postage- 
stamp images, both of which are integral to later discus- 
sions. There is no detectable counterpart at the X-ray 
position of CXOU J141312.3-652013 (hereafter CG X-1). 
CXOU J141310.0-652044 (hereafter CG X-2) coincides with 
an my = 22.2 companion embedded in a probable H II re- 
gion, as indicated by the diffuse emission surrounding the op- 
tical point source in the HST Ha image in Figure 2. The 
bright, compact Ha emission (wea = 17.9) from this source 
may originate in either a compact H II region or a young SNR. 
In general, we found that many of the serendipitous X-ray 
sources appear to be coincident with optical star-forming re- 
gions: CXOU J141310.3-652017 with the inner 40 pc star- 
forming ring; CXOU J141310.1-652029 with the outer 200 pc 
star-forming ring; CG X-2 with an inter-arm Ha complex; and 
CXOU J141309.2-652017, CXOU J141310.4-652022, CXOU 
J141312.2-652007 and CG X-1 with the spiral arms. 

After assessing the optical properties of the X-ray point 
sources, we used X-ray-to-optical flux ratios to evaluate dif- 
ferent X-ray emission mechanisms. For instance, Mac- 
cacaro et al. (1988) found that normal stars have -4.3 < 
\og(Fx/F\j) < -0.3 while background galaxies and AGN have 
-1.5 < log(Fx/Fy) < 1-9.''' RS CVn stars, XBs and SNRs are 
harder to characterize because their logiFx /Fy) ratios can range 
from values close to those of normal stars (RS CVn stars) up to 
ratios greater than 3.0 (XBs; see Maccacaro et al. 1988). This 
large range in X-ray-to-optical flux ratio leaves some ambiguity 
in the classification scheme for sources with log(Fx/Fy) < 1.9, 
but the scheme is fairly clear cut for sources with larger ratios. 
In the case of Circinus, the off-nuclear sources have log(Fx/Fv) 
ranging from -0.9 to 3.4 (see column 7 of Table 1). We find that 
five off-nuclear point sources are unambigiously classified as 
XBs or SNRs. The remaining sources could be similarly clas- 
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sified, but the X-ray-to-optical flux ratios alone cannot exclude 
the possibility that eight sources could be background galaxies 
or AGN and two sources could be normal stars. We address the 
probability of finding foreground or background sources coin- 
cident with Circinus below. 

Nearly all of the detected X-ray point sources are coinci- 
dent with the disk of Circinus and lie within 2' of the nu- 
cleus. However, Circinus lies behind a crowded stellar re- 
gion at low Galactic latitude (b = -3.° 8), so the chance super- 
position of Galactic X-ray sources is a vaUd concern. Fur- 
thermore, Circinus has a Galactic longitude of / = 311.°3, so 
our line-of-sight to Circinus intersects a sizable fraction of the 
Galaxy. To estimate the number of foreground or background 
sources expected within the central region of Circinus, we used 
a deep Chandra observation of the Galactic Ridge (/ = 30°, 
b = 0°; K. Ebisawa et al., in preparation). This observation 
is at a similar Galactic latitude and has a limiting X-ray flux of 
^2-10 kev = 3 X 10"''' erg cm"^ s"' . Since Circinus and the Galac- 
tic Ridge have comparable angular distances from the Galactic 
Center, they should also have comparable source-count distri- 
butions. In the 2-10 keV band, the Circinus zeroth-order image 
has a flux threshold of F2-10 keV = 6 x 10"'^ erg cm"^ s"' assum- 
ing r = 1 .76 and A^h = 1 -5 x 10^^ cm"^. From the Galactic Ridge 
measurements, we expect to detect « 300 sources per deg^ in 
a 67 ks Chandra HETGS observation, corresponding to « 1.0 
foreground or background sources within a 2' radius of Circi- 
nus at our detection threshold (see Figure 1). Therefore, most 
of the sources within 2' of the nucleus should be associated with 
Circinus. 

Alternatively, we could ask how many sources we would ex- 
pect to find in the field outside the disk of Circinus (i.e., > 
2' from the Circinus nucleus)? If we exclude Circinus from 
the ACIS S3 chip subarray, we expect to find « 2 sources 
in the remaining area (w 27 arcmin^) after correcting for vi- 
gnetting losses. Only one source is detected within this re- 
gion (CXOU J141257. 8-651425), indicating that our use of 
the Galactic Ridge number counts is appropriate and perhaps 
somewhat pessimistic. 

Finally, we note that sources significantly brighter than our 
limiting flux have much more stringent limits. Since the Chan- 
dra Galactic Ridge number counts are poorly sampled above 
■P2-10 keV = 1 X 10"'^ erg cm"^ s"', we combined them with the 
results of thsASCA Galactic Plane Survey (Sugizaki et al. 1999) 
to estimate that there is a < 0.06% chance that either of the 
bright off-nuclear sources shown in Figure 2 (i.e., CG X-1 and 
CG X-2) are foreground or background X-ray sources. 

4. TIMING ANALYSIS 

4.1. Short- Term Variability 

The Chandra observation of Circinus was of sufficient du- 
ration to evaluate short-term timing characteristics for many 
of the off-nuclear sources. To determine objectively the exis- 
tence of significant variations in the count rate, we used the 
Kolmogorov-Smirnov (KS) statistic to test the null hypothe- 
sis that each source plus background rate was constant over 
the duration of the exposure.'^ Thus a source with a low 
KS probability (Peon) has a high probability of being variable 
(1 - Peon)- Since the timing gap in the observation could lead 
to inaccurate variability estimates, we modified the KS test 
to exclude the 1.3 ks gap between the continuous 6.9 ks and 
60.2 ks segments of the observation. We find that sources 

" The data were not binned for the KS test. The minimum time resolution was the 
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Fig. 3. — Light curve of CG X-1 over the entire duration of the observation 
(i.e., 6.9 ks + 1.3 ks gap + 60.2 ks) binned in 2 ks intervals. The period is 
27.0 ± 0.7 ks. In that time, the flux varies by more than a factor of 20. The 
1.3 ks gap (shaded band) was interpolated using the count rates found from av- 
eraging the 1 ks data segments before and after the gap. The eiTors bars have 
been calculated following the method of Gehrels (1986). 

CXOU J141312.2-652007, CG X-1, and CXOU 
J141312.6-652052 varied at the > 99% confidence 
level, while sources CXOU J141307.5-652106 and 
CXOU J141310. 4-652022 were less convincing, varying at 
the K, 92% and « 98% confidence levels, respectively. Al- 
though CXOU J141307.5-652106 only has a KS probabiHty of 
92% and 7.7 background-subtracted counts, all of these counts 
were detected in a contiguous 20 ks segment of the observa- 
tion, indicating it is a strong flaring source. In total, three of 
these short-term variable sources appeared to flare on ^ 20 ks 
timescales, while one exhibited a steady increase in flux over 
the course of the observation and another displayed periodic 
behavior (see below). The short-term KS statistic results for all 
16 Chandra sources and the nature of variability observed are 
summarized in columns 3 and 10 of Table 1, respectively. 

Many of our serendipitous X-ray sources are likely to be XBs 
(see §3). Since a significant fraction of XBs are X-ray pulsars 
with pulse periods ranging from 0.01-1,000 s (e.g.. White, Na- 
gase, & Parmar 1995), we generated power spectra for all of the 
Circinus X-ray sources to search for pulsations. These power 
spectra also allowed us to search for other periodic phenomena 
(e.g., eclipses or regular X-ray bursts). Only one source, CG X- 
1, was found to display periodic behavior above the frame-time 
limit of 2.1 s. Figure 3 presents the light curve of CG X-1, 
which shows periodic flux variability by at least a factor of 20. 
Because the 1.3 ks gap could introduce a small artificial dip 
in the light curve, we interpolated the flux values during this 
time using the count rates found from the 1 ks data segments 
before and after the gap. We used the period-fitting routine of 
Stetson (1996) to quantify the period and its one-sigma error as 
7.5 ±0.2 hr. To study the light curve in more detail, we folded 
it about its period as shown in Figure 4 to increase the signal- 
to-noise. The flux of the source exhibits a gradual rise early in 
its phase (~ 20% of cycle), then flattens out near the peak flux 
40% of cycle), and abruptly fades to a minimum near zero 
(^ 40% of cycle) to close out the period. As the flux gradually 
increases in Figure 4 the spectrum hardens briefly, suggesting 
that the emission during this portion of the period may be par- 
tially absorbed at soft energies. 

The simplest explanation for the variability of CG X-1 is 
an eclipse by an evolving binary companion, although a pure 

frame time. 
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Fig. 4.— (Top) Light curve of CG X-1 folded about a period of 27.0 ks 
and binned in 1 ks intervals. The error bars have been calculated following the 
method of Gehrels (1986) assuming only photon noise. (Bottom) Hardness 
ratio (HR) of CG X-1 folded about a period of 27.0 ks and binned in 3 ks in- 
tervals. HR = (H - S)/(H + S) where H = 2-10 keV counts and S = 0.5-2 keV 
counts. The errors bars for this quantity have been calculated following the 
"numerical method" described in §1.7.3 of Lyons (1991). Prior to folding the 
flux and HR light curves, the 1.3 ks gap was interpolated using the count rates 
found from averaging the 1 ks data segments before and after the gap. 

eclipse model has difficulty explaining certain aspects of the 
light curve. For instance, the fraction of the period spent 
in eclipse is related to the size of the Roche Lobe of the 
binary companion and hence to the companion-to-compact- 
object mass ratio. We find that the eclipse fraction of CG X-1 
is w 40%, implying a mass ratio of w 30 (see the Appendix 
of Pringle & Wade 1985). Such a large mass ratio is diffi- 
cult to reconcile with other observed quantities. For instance, 
if CG X-1 resides in the disk of Circinus at a distance of 3.8 
Mpc (see §3), then its average and peak X-ray luminosities are 
3.7 X 10^^ erg s~' and 6.6 x 10^^ erg s"'. The peak luminosity 
implies a compact-object mass of M > 50Mq if it is radiat- 
ing at or below its Eddington luminosity. The mass limit for 
the compact object and the mass ratio predicted by a simple 
eclipse model are incompatible with the standard model of stel- 
lar evolution (i.e., a companion mass of > 15OOM0 would be 
required). Therefore, in the context of an eclipse model, the 
large eclipse fraction implies that this binary system must ei- 
ther be experiencing strong dynamical Roche Lobe overflow, 
or, alternatively, the X-ray emission region must be at least par- 
tially occulted by some other means such as an accretion disk 
around the compact object. 

Another possible scenario for generating the large-amplitude 
periodic variations seen from CG X-1 is by modulation of the 
accretion rate (e.g., from accretion-disk instabilities). One such 
example is the extremely variable and luminous Galactic X-ray 
source, GRS 1915+105. In its active state, the observed X-ray 
intensity of GRS 1915+105 is found to vary dramatically on 
timescales ranging from seconds to days, occasionally display- 
ing quasi-periodic episodes (e.g., Greiner, Morgan, & Remil- 
lard 1996). Some of these quasi-periodic episodes coarsely re- 



semble the light curve of CG X-1 (although the timescales are 
generally shorter for GRS 1915+105). Thus CG X-1 could also 
plausibly be an analog of GRS 1915+105 within Circinus. 

We note that the asymmetric rise and decline in intensity of 
the CG X-1 light curve resembles those of AM Her type sys- 
tems (e.g., Warner 1995). This class of CV typically has pe- 
riods of 1-4 hr, intrinsic X-ray luminosities between lO-'' < 
Lq.s-iq keV < 10"'"' erg s"', and binary companions occupying a 
narrow range of spectral types (M2V-M6V; Cordova 1995; Ver- 
bunt & van den Heuvel 1995; Warner 1995). These characteris- 
tics restrict the X-ray-to-optical flux ratio to 3 < Fx/Fy < 300. 
Given the X-ray flux associated with CG X-1 (Table 1), we 
would expect an M star companion to have my « 16-20, which 
is easily ruled out by the //ST limit. The //ST limit requires that 
any M2V-M6V companion star of CG X-1 be at least 1.2 kpc 
away. The implied X-ray luminosity (> lO"*^ erg s"') at this dis- 
tance is at least 10 times larger than the most luminous known 
AM Her system. Furthermore, the period (7.5 hr) is nearly two 
times longer than the longest period AM Her system known, 
and the X-ray spectrum is inconsistent with the two-component 
model typical of these systems (see §5). Based on these facts 
alone, a Galactic identification is doubtful. Moreover, the ob- 
served X-ray flux of this object limits the probability of it being 
a foreground source to < 0.06% (see §3). We also note that 
AM Her systems in particular comprise only a small fraction 
of all Galactic X-ray emitters (-- 2%; e.g., Motch et al. 1998). 
Therefore the most promising explanantion for CG X-1 is an 
intermediate-mass black hole within Circinus. 

4.2. Long-Term Variability 

Circinus was also observed with the ROSAT HRI (« 5" 
FWHM; 0.1-2.4 keV) in 1995 September (Guainazzi et al. 
1999). With a total exposure time of only 4.2 ks, this 
HRI observation is limited to source detections above 1.7 x 
10""* erg cm-2 s"' (0.5-2.0 keV). CG X-1 was detected with 
a 0.5-2.0 keV flux of (8.7 ± 2.9) x 10"'* erg cm'^ s"', a fac- 
tor of 2 lower than the average 0.5-2.0 keV flux mea- 
sured with Chandra (1.7 x 10"^^ erg cm"^ s~'), but still well 
within the range of intensity variations observed by Chandra. 
CG X-2, on the other hand, is completely absent from the 
HRI image, corresponding to a 95% confidence upper limit 
of 3.3 X lO"'** erg cm"^ s"' (using the method of Kraft, Bur- 
rows, & Nousek 1991) in the 0.5-2.0 keV band. This is 
a factor of at least 11 smaller than the 0.5-2.0 keV flux of 
(4. 1 ±0. 1) X 10"' erg cm"^ s"' measured with Chandra. CG X- 
2 is coincident with an H II region, and the Chandra radial 
profile is consistent with that of a point source, indicating that 
this source could be either a moderately young SNR or tran- 
sient XB in Circinus. Observations taken with the Australian 
Telescope Compact Array (ATCA) in 1995 with a beam size of 
^Beam = 0."9 X 0."8 show an unresolved radio source (S^ cm ^ 
1.2 mJy, aj2ooo = Mi^n^lO^O (5j2ooo = -65°20'44"; M. Elmout- 
tie, private communication) coincident to within < 0"6 of the 
X-ray source, while none of the other radio point sources coin- 
cide with X-ray sources (see Figure 8 of Elmouttie et al. 1998). 
Given that there are only a few significant radio point sources 
in the Circinus field at a wavelength of 3 cm, the probability of 
this pair being a random alignment is < 0. 1 %. The combination 
of bright, point-like radio, optical, and X-ray emission strongly 
suggests that this object is a young SNR or an H II region. We 
can rule out radio emission from a bright H II region, however, 
since the observed flux at 3 cm is approximately two orders of 
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magnitude higher than expected."^ A young SNR classifica- 
tion for CG X-2 is further supported by our spectral fitting in 
§5.1. The detection of a moderately strong radio source in 1 995 
and the lack of a bright X-ray source in 1996 is curious but not 
implausible given that the X-ray and radio emission are likely 
to result from physically distinct and separate emission regions 
within a young SNR (e.g., Schlegel 1995). Ongoing optical 
monitoring programs of Circinus over the past eight years have 
failed to detect any bright supernovae down to my « 15 (R. 
Evans 2000, private communication; A. Williams 2000, private 
communication). Thus, if CG X-2 exploded after 1993, it must 
have been optically faint. 

5. SPECTRA 
5.1. Zeroth-order 

The ACIS zeroth-order spectra of the off-nuclear sources 
were extracted for sources with more than w 90 counts using 
either the 90% or 95% encircled-energy extraction radii as de- 
scribed in §3. Circinus was placed near the node 0/1 boundary 
on the S3 chip, so several sources including the nucleus were 
dithered between nodes and 1 . Since the variation in the en- 
ergy response of the S3 chip across these nodes is large, we 
chose to extract the spectra and response matrices for each node 
separately and fit both spectra with one model in XSPEC. To as- 
sess the spectral nature of these sources, we initially fitted the 
spectra with absorbed power-law models. We found that only 
one off-nuclear source required more extensive spectral model- 
ing (CG X-2). The best-fit models with 90% confidence errors 
(Ax^ = 2.7) are listed in columns 8 and 9 of Table 1. 

The range of spectral properties for these objects is quite di- 
verse (e.g., r ~ 0.7-2.7). All of these off-nuclear sources ex- 
hibit spectral cutoffs below ~ 1 keV that are best fitted with 
column densities larger than the Galactic value (see §1). These 
large absorption values are consistent with the fact that many of 
the off-nuclear X-ray sources appear to reside in star-forming 
regions and spiral arms in Circinus (see §3). The spectra of the 
two brightest off-nuclear sources have adequate counts to con- 
strain physical models of emission and are discussed in more 
detail below. 

In Figure 5, we present the spectrum of the periodic source 
CG X-1. CG X-1 happened to fall near the boundary of 
node and node 1, such that w 70% of the photons from 
this source were dithered onto node and « 30% onto node 
1. We show the node spectrum in Figure 5. This spec- 
trum was initially fitted with a moderately soft power law 
with r = 2.40!|]:22 and A'h = il-26^ jl) x lO^^ cm-^ (x'^ = 75.4 
for 59 degrees of freedom). Although this model proved ac- 
ceptable, Makishima et al. (2000) have shown that the X-ray 
spectra of super-Eddington sources in spiral galaxies are often 
successfully described by multicolor disk blackbody emission 
(MCD; Mitsuda et al. 1984) arising from optically thick stan- 
dard accretion disks around black holes. We found that a MCD 
model with an inner-disk temperature of Jin = 1 •35;|;q Jj keV and 
A^H = (6.39![;^^) X 10^^ cm"^ worked somewhat better than the 
power-law model (x^ = 61.0 for 59 degrees of freedom) and 
that the derived spectral parameters agree well with typical in- 
nermost disk temperatures derived for other super-Eddington 
sources (i.e., = 1.1-1.8 keV; Makishima et al. 2000). 
Thus, we find that the spectral properties of this source are 

We assume that a typical bright H II region has an electron temperature of 10'* K, 
see Osterbrock 1989). 

" In most cases, a Raymond-Smith thermal plasma model (Raymond & Smith 1977) 
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Fig. 5. — (Top) The upper panel shows the ACIS-S node spectrum of 
CG X-1, modeled with a MCD model with an inner-disk radius temperature of 
Tin = 1.35+55 13 keVand/VH = (6.39+||5) X lO" cm'^. The lower panel shows 
the residuals of the fit measured using the statistic. (Bottom) Confidence 
contours between and Nh for CG X-1. The contours represent the 68%, 
90% and 99% confidence levels for two parameters of interest. 

fully consistent with X-ray emission from an intermediate-mass 
black hole candidate in the disk of Circinus. We also attempted 
to fit the spectrum of CG X-1 assuming the typical spectral 
model for an AM Her system (i.e., the most likely Galactic 
candidate) to remove any remaining concerns about a possi- 
ble Galactic nature of CG X-1. Spectra of AM Her systems 
have substantial 6.4 keV iron Ka emission lines and continua 
comprised of a blackbody component with kT^B ~ 20 eV and a 
bremsstrahlung component ^7br ~ 30 keV (see §6.4 of Warner 
1995). Such a model does not fit the observed spectrum well 
and is rejected with ^ 99% confidence. 

The other bright X-ray source, CG X-2, was initially fitted 
with a power law with T = 1.77;|^ 22 and an absorption col- 
umn of A^H = X 10^^ cm-2 (see Figure 6). How- 
ever, this fit was statistically unacceptable with = 106.5 
for 74 degrees of freedom and a null-hypothesis probabil- 
ity of Pix^lu) = 8.2 X 10"^. The spectrum shows large pos- 
itive residuals around 6.9 keV. The addition of a Gaussian 
emission line at 6.89;!;qq9 keV with an intrinsic line width of 
190!^^ eV improved the fit significantly (x^ = 63.2 for 71 

electron and proton densities of 10^ cm~^, a size of less than 10 pc (for details 

was equally acceptable. 
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Fig. 6. — (Top) The upper panel shows the ACIS-S node spectrum of 
CG X-2, modeled with a neutral hydrogen column of A^h = (l-lS^gl^) x 10^2 
cm~2 and a power law with F = 1 -77;^^ Ij . The lower panel shows the residuals 
of the fit measured using the statistic. Note the large residuals around 
6.67-6.97 keV. (Bottom) The confidence contours for fit parameters of a single 
Gaussian emission line at 6.89 keV. The contours represent the 68%, 90% and 
99% confidence levels for two parameters of interest. 

degrees of freedom). The measured equivalent width of the line 
is 1 .59;!;q 4g keV. We note that this source is offset from the nu- 
cleus by 23" and totally dominates over the faint background. 
Furthermore, the energy of this line (6.89 keV) is quite different 
from the emission-line feature found in the rest of the reflection 
halo (« 6.4 keV, for details see Sambruna et al. 2001a) and is 
not associated with that component. 

To confirm that the central energy and line width of the 
6.89 keV line were not due to possible calibration problems, 
we measured the line energies and widths of the 5.89 keV and 
6.49 keV caUbration Unes in a caUbration dataset closest to our 
observation date (obs62017, 2000 June 18, -120° C). From the 
calibration data, we extracted a large region on node of the 
S3 chip near the position of CG X-2 and fit two Gaussian emis- 
sion Unes to the calibration lines using the same response ma- 
trix as used for CG X-2. The best-fit values were line energies 
of 5.94!U keV and 6.54!g;J]j keV and line widths of 53!|]^ eV 
and 56;|;q3 eV, respectively. Thus there appears to be a shift of 
~ 50 eV in the central fine energy and uncertainty in the line 
width of « 60 eV. The true central energy of the 6.89 keV line 
should therefore be 6.84 keV. This 6.84 keV line is also broader 
than the instrumental resolution of ACIS, suggesting that it is 
either intrinsically broad or a blend of emission from two or 



more lines. 

Given the central energy of the 6.84 keV emission line, it 
is likely to be a blend of iron K-sheU transitions at 6.67 keV 
(He-like line complex) and 6.97 keV (H-like line). To test this 
hypothesis, we modeled the spectrum of CG X-2 using two 
Gaussian emission components and an absorbed power law. We 
fixed the ratios of the He-like complex and H-like line energies 
and fine widths to be 0.957 and 1.0, respectively, but we al- 
lowed the energy and width of the combination to vary. The 
normalizations of the two lines were allowed to vary indepen- 
dently. This model improved the to 57.5 for 70 degrees of 
freedom and is better than the single Gaussian model at » 99% 
confidence using the F-test. The best-fit energies of the iron 
lines are 6.74^q "3 keV and 7.05;|;o q3 keV and are consistent with 
the line rest energies after removing the 50 eV blueshift due to 
calibration errors (see above). The intrinsic fine width of each 
line is 70;!;27 eV and is consistent with unbroadened emission 
lines (compare to line widths of calibration lines above). Al- 
lowing the energies and widths of the fines to vary freely does 
not significantly change the best-fit results, although the errors 
associated with each parameter are much larger. 

There are two scenarios that could explain such large equiva- 
lent width iron K-shell line emission. We first explore the possi- 
bihty that this source is an X-ray binary where the strong emis- 
sion lines must result from fluorescence and resonant scattering 
either in a wind or a corona around the accretion disk (e.g., see 
§1.6.2.3 of White, Nagase, & Parmar 1995). In this scenario, 
the direct line-of-sight to the X-ray source must be blocked; 
otherwise the continuum emission from the source would dilute 
the line emission and reduce its equivalent width below what 
is observed. Scattering of the continuum by electrons is rela- 
tively inefficient (« 1-5%), so the intrinsic X-ray luminosity of 
this source would have to be a factor of « 20-100 times larger 
than observed. We also detect no short-term variabihty for this 
source, so the direct continuum would have to be blocked for 
more than 68 ks. Since the luminosity of this source is already 
^0.5-10.0 keV = 3.4 X lO'^'^ erg s"^ and exhibits no short-term vari- 
ability, we consider an X-ray binary scenario for CG X-2 to be 
unlikely. 

The second scenario to explain the X-ray line emission is that 
of a young SNR (see §4.2). Within this scenario, the line emis- 
sion could be generated from either the initial forward shock 
that passes through the outer layers of the star and propagates 
into the circumstellar material or from the high-emissivity re- 
verse shock that develops at later times. The large X-ray lu- 
minosity of CG X-2 is comparable to those of other young (< 
20 yr old) SNRs detected at X-ray wavelengths (e.g., Schlegel 
1995; Immler, Pietsch, & Aschenbach 1998a, b; Fox et al. 2000; 
Pooley et al. 2001), adding merit to such a classification. To fur- 
ther test this hypothesis, we fitted the spectrum with Raymond- 
Smith thermal plasma and non-equilibrium ionization colli- 
sional plasma models (Borkowski, Sarazin, & Blondin 1994). 
Figure 7 shows the best-fit model for an absorbed Raymond- 
Smith thermal plasma [kT = 10.82!|^^ keV, Z = (1.35!°j4)ZQ, 
A^H = (9.21+};^^) X 10^^ cm"2]. While this model provides an ac- 
ceptable fit to the data (x^ = 77.1 for 73 degrees of freedom), it 
clearly underpredicts the width of the highly ionized iron emis- 
sion, leaving a gap between the H-like and He-like iron lines 
around 6.8 keV. The non-equilibrium ionization model proved 
no better at fitting the broad iron-line blend. 

Only w 12 young SNRs have X-ray detections at all, 
and when adequate ROSAT, ASCA or Chandra spectra are 
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Fig. 7. — The upper panel shows the ACIS-S node spectrum of CG X- 
2, modeled with a neutral hydrogen column of A'h = (9.21t} 4g) X 10^^ cm"^, 
a Raymond-Smith thermal plasma model with kT = 10.82|';2 12 keV and Z = 
(1.35;*^ 54)Zq. The lower panel shows the residuals of the fit measured using 
the statistic. The emission-line blend centered at 6.89 keV is poorly fit. 

available, most display high thermal temperatures {kT = 5- 
10 keV) similar to CG X-2. Furthermore, only five young SNRs 
have been observed with adequate sensitivity above 5 keV 
to detect any K-shell iron emission. Three of these sources 
(SN 1986J in NGC 891, Houck et al. 1998; SN 1993J in M81, 
Kohmura et al. 1994; SN 1998S in NGC 3877, Pooley et al. 
2001) appear to have substantial 6.8 keV iron-line emission 
with equivalent widths of « 0.3-2 keV. Thus the spectrum of 
CG X-2 appears to be consistent with the X-ray emission seen 
from other young SNRs. The failure of simple physical models 
(e.g., a Raymond-Smith model) to fully describe the iron K- 
shell emission from CG X-2 is intriguing, but not implausible 
considering the scarcity of young, nearby SNRs and the lack 
of detailed X-ray spectra. In addition to highly ionized iron K- 
shell emission, we might expect to detect faint line emission 
from other highly ionized elements such as sulphur and silicon. 
The spectrum of CG X-2, however, is too noisy to place reliable 
limits on such emission. 

5.2. HEG/MEG 

For the two brightest off-nuclear sources, CG X-1 and CG X- 
2, we extracted MEG and HEG grating spectra (see Sambruna 
et al. 2001b for details of the extraction procedure). We found 
that the high-resolution spectra were consistent with the ACIS- 
S spectra, but that the statistics were generally poor (less than 
3(7 line detections in the case of CG X-2). 

6. LUMINOSITY FUNCTION OF RESOLVED SOURCES 

Luminosity functions (LFs) are often a useful tool for shed- 
ding physical insight into the properties of source populations. 
Figure 8 shows the cumulative distribution of the Circinus off- 
nuclear sources as a function of luminosity. For comparison, 
we also show three other galaxies spanning a large range of 
star formation rates for which similar 0.5-10.0 keV source 
distributions could be obtained: the Small Magellanic Cloud 
(SMC; Yokogawa et al. 2000), M82 (Griffiths et al. 2000) and 
NGC 3256 (R Lira et al., in preparation). The LFs of these 
three galaxies all appear to follow the same trend such that 

oc E^'^^ ■ However, the LF of Circinus shows a distinct "kink" 
around 10^^ - 10^^ erg s"' that is not evident in the other LFs. 
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Fig. 8. — Cumulative number of sources versus 0.5-10.0 keV luminosity for 
Circinus (solid), the SMC (dotted), M82 (short-dashed), and NGC 3256 (long- 
dashed). The sources in the SMC, M82 and NGC 3256 all appear to increase 
with decreasing luminosity such that N oc L^'^^ ■ The Circinus sources do not 
follow this trend (see §6). 

One possible explanation for this anomalous LF may be the 
large attenuation of soft X-ray photons from both the large 
Galactic column and absorption in the HETGS. For example, 
if the Circinus source population has a LF slope similar to that 
of the other galaxies, but the low-luminosity sources in Circi- 
nus are substantially softer than the high-luminosity ones (an 
effect seen in NGC 4038/4039; A. Zezas et al. 2001, in prepa- 
ration), a large fraction of the faint, soft X-ray sources might go 
undetected. 

To test this hypothesis, we assumed that, in the absence of 
Galactic and HETGS absorption, the LF of Circinus has the 
same average intrinsic slope as that of the other galaxies in Fig- 
ure 8 (i.e., oc Lx"'^^)- 

Since the brightest source detected in 
Circinus has a 0.5-10.0 keV luminosity of 3.4 x lO-'^ erg s"', 
we should expect to detect 28 unabsorbed sources at the lu- 
minosity limit of our HETGS observation (i.e., Lo.s-io.okev = 
2 X lO"'^ erg s"'). For simplicity, we assumed that all sources 
below L().5_io,() keV = 1 X 10^^** erg s"' were well-fit by a power- 
law spectrum with an average photon index F. We then used 
PIMMS (Mukai 2000) to estimate the number of sources that 
would fall below our flux threshold as we varied the average F 
from 1.76 (i.e., what we find for the bright sources) to softer 
values, assuming an HETGS observation and an average neu- 
tral hydrogen absorption column of A^h = 1.5 x 10^^ cm"^. We 
found that in order to reproduce the observed 15 off-nuclear X- 
ray sources, an average photon index of F w 3.2 is necessary 
for the fainter sources. 

While this hypothesis provides a viable explanation for the 
"kink" seen in the LF of Circinus, the average spectral index 
required for the fainter sources is significantly softer than the 
softest spectral index directly measured via spectral fitting anal- 
ysis of the brightest sources in the field (see § 5.1 and Table 1). 
This large contrast argues against such a hypothesis, but we 
note that a similar trend of systematic softening, although con- 
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siderably milder, is seen in Chandra point-source population 
studies of nearby galaxies such as M81 (Tennant et al. 2001), 
NGC 4038/4039 (A. Zezas et al. 2001, in preparation), and 
NGC 4647 (Sarazin, Irwin, & Bregman 2000). The small num- 
ber of counts in these studies, however, allow only crude hard- 
ness ratio determinations. Such a trend is not found in a deep 
XMM-Newton study of M31 (Shirey et al. 2001), where there 
are adequate counts to perform detailed spectral analysis on all 
sources above io.s-io.o kev « 1 x 10^^ erg s~^ Unfortunately, 
a large statistical sample is not yet available to clarify whether 
this trend of systematic spectral softening is common among 
the fainter point sources in nearby galaxies or whether an alter- 
nate explanation for LF of Circinus exists (e.g., there is some in- 
trinsic difference between Circinus and the other galaxies such 
as morphological type, gas mass, star formation rate, star for- 
mation history, etc. which could explain the "kink"). 

7. DISCUSSION AND SUMMARY 

We have performed a systematic analysis of 16 point sources 
coincident with the disk of Circinus in a Chandra HETGS 
observation, down to a limiting X-ray luminosity of 2 x 
10^^ erg s"'. We find that 

• Nearly half of the sources coincide with the 
star-forming features of the galaxy, and only four 
off-nuclear sources have optical counterparts. 

• Four sources exhibit short-term (~ 1 day) variability, 
and one source exhibits only long-term (~ 4 yr) 
variability. 

• The X-ray spectra of these objects are best-fit by a large 
absorption column above the Galactic value and a range 
of power-law photon indices from F = 0.7-2.7. 

Based on the large observed absorption columns, spatial coinci- 
dence with the Circinus galaxy, and large X-ray-to-optical flux 
ratios, we conclude that the bulk of these off-nuclear sources 
are XBs and/or luminous SNRs located within Circinus. The 
shape of the Circinus LF appears to be different from those of 
other galaxies and may indicate a bias against faint soft X-ray 
sources introduced by the large Galactic column and absorption 
in the HETGS. Any definite conclusions about the Circinus LF, 
however, should wait until a large, complete sample of nearby 
galaxies can be assembled for comparision. 

The two most luminous off-nuclear sources, CG X-1 and 
CG X-2, are remarkable. The distinct periodic variability, high 
X-ray luminosity, and soft X-ray spectrum of CG X-1 provide 
strong support for an interpretation as a > 50 black hole 
in an accreting binary system in Circinus. High-luminosity X- 
ray sources such as this one are fairly common among nearby 



galaxies (see §1), but the physical nature of these objects is of- 
ten poorly determined because of the limited photon statistics. 
The X-ray properties of the periodic source CG X-1 provide 
the strongest evidence to date that at least some of the super- 
Eddington X-ray sources detected in other galaxies are indi- 
vidual XBs with large inferred black-hole masses, challenging 
both stellar evolution and accretion disk models (e.g., Taniguchi 
et al. 2000; Watarai, Mizuno, & Mineshige 2001). Unfortu- 
nately, the lack of an optically detectable stellar companion 
limits the degree to which we can constrain the physical pa- 
rameters of this system. Further X-ray observations of CG X-1 
are important to tighten constraints on the variabiUty of its pe- 
riod, luminosity and spectrum. If this source indeed undergoes 
periodic eclipses, then the period and shape of the light curve 
should remain fairly stable. However, if the apparent periodic- 
ity is due to accretion disk instabilities, one might expect the 
X-ray hght curve to "evolve" away from what we see in Fig- 
ure 3. 

Likewise, the large X-ray luminosity, the presence of a very 
large equivalent width iron K-shell emission line, the hard X- 
ray spectrum, and the association with both a radio point source 
and an Ha-detected H n region make CG X-2 a prime candidate 
for a young SNR in Circinus. Continued X-ray monitoring of 
this source will determine how its luminosity and spectral prop- 
erties evolve with time. Nearby, X-ray-bright supemovae like 
CG X-2 are relatively rare (Schlegel 1995; Houck et al. 1998), 
but they represent an excellent opportunity to study the physical 
interaction of a supernova with a dense circumstellar environ- 
ment. 

The ubiquity of super-Eddington off-nuclear X-ray sources 
in nearby galaxies has become an established but poorly un- 
derstood phenomenon. Our findings for Circinus support the 
hypothesis that the off-nuclear, super-Eddington X-ray sources 
in nearby galaxies are the product of both young SNRs and 
intermediate-mass black holes in XBs. 
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